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We explore the feasibility and astrophysical consequences of a new long-range (7(1) gauge field 
("dark electromagnetism" ) that couples only to dark matter, not to the Standard Model. The 
dark matter consists of an equal number of positive and negative charges under the new force, 
but annihilations are suppressed if the dark matter mass is sufficiently high and the dark fine- 
structure constant a is sufficiently small. The correct relic abundance can be obtained if the dark 
matter also couples to the conventional weak interactions, and we verify that this is consistent 
with particle-physics constraints. The primary limit on a comes from the demand that the dark 
matter be effectively collisionless in galactic dynamics, which implies a < 10 -3 for TeV-scale dark 
matter. These values are easily compatible with constraints from structure formation and primordial 
nucleosynthesis. We raise the prospect of interesting new plasma effects in dark matter dynamics, 
which remain to be explored. 

I. INTRODUCTION 

A wide variety of cosmological observations seem to point to a two-component dark sector, in which approximately 
73% of the energy density of the universe is in dark energy and 23% is in non-baryonic dark matter (DM) . Ordinary 
matter constitutes the remaining 4% [1 . The physics of the dark matter sector is plausibly quite minimal: an excellent 
fit to the data is obtained by assuming that dark matter is a cold, collisionless relic, with only the relic abundance as 
a free parameter. The well-known "WIMP miracle" 0[3J[1] is the fact that a stable, neutral particle with weak-scale 
mass and coupling naturally provides a reasonable energy density in DM. Particles of this type arise in models of low- 
scale supersymmetry or large extra dimensions [5J, and provide compelling DM candidates. In the contemporary 
universe, they would be collisionless as far as any conceivable dynamical effects are concerned. 

Nevertheless, it is also possible to imagine a rich phenomenology within the dark sector. The dark matter could be 
coupled to a relatively strong short-range force that could have interesting consequences for structure on small scales 
[HIE]. Alternatively, DM could also be weakly coupled to long-range forces, which might be related to dark energy 
[5]. One difficulty with the latter is that such forces are typically mediated by scalar fields, and it is very hard to 
construct natural models in which the scalar field remains massless (to provide a long-range force) while interacting 
with the DM at an interesting strength. 

In this paper, we explore the possibility of a long-range gauge force coupled to DM, in the form of a new unbroken 
abelian field, dubbed the J7(1)d "dark photon." We imagine that this new gauge boson 7 couples to a DM fermion 
X, but not directly to any Standard Model (SM) fields. Our model is effectively parameterized by only two numbers: 
m x , the mass of the DM, and d, the dark fine-structure constant. If m x is sufficiently large and a is sufficiently small, 
annihilations of DM particles through the new force freeze out in the early universe and are negligible today, despite 
there being equal numbers of positively- and negatively-charged particles. The dark matter in our model is therefore 
a plasma, which could conceivably lead to interesting collective effects in the DM dynamics. 

Remarkably, the allowed values of m x and a seem quite reasonable. We find that the most relevant constraint 
comes from demanding that accumulated soft scatterings do not appreciably perturb the motion of DM particles in 
a galaxy over the lifetime of the universe, which can be satisfied by a ~ 10~ 3 and m x ~ TeV. For values near these 
bounds, the alterations in DM halo shapes may in fact lead to closer agreement with observation [6]. However, for such 
regions of parameter space, if U (1) d were the only interaction felt by the \ particles, the resulting relic abundances 
would be slightly too large, so we need to invoke an additional annihilation channel. We show that \ can m f& c t be 
a WIMP, possessing SU(2)l quantum numbers in addition to U(1)d charge. Such a model provides the correct relic 
abundance, and is consistent with particle-physics constraints so long as the mixing between ordinary photons and 
dark photons is sufficiently small. 

We consider a number of other possible observational limits on dark electromagnetism, and show that they do not 
appreciably constrain the parameter space. Since the DM halo is overall neutral under U(1)d, there is no net long- 
range force that violates the equivalence principle. Although there are new light degrees of freedom, their temperature 
is naturally lower than that of the SM plasma, thereby avoiding constraints from Big-Bang Nucleosynthesis (BBN) . 
Energy loss through dark bremsstrahlung radiation is less important than the soft-scattering effects already mentioned. 
The coupling of DM to the dark radiation background can in principle suppress the growth of structure on small scales, 
but we show that the DM decouples from the dark radiation at an extremely high redshift. On the other hand, we 
find that there are plasma instabilities {e.g. the Weibel instability) that can potentially play an important role in the 
assembly of galactic halos; however, a detailed analysis of these effects is beyond the scope of this work. 
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The idea of an unbroken U(l) coupled to dark matter is not new. 1 De Rujula et al. [TO] explored the possibility 
that dark matter was charged under conventional electromagnetism (see also [II] Q2J [TO] [14] ) ■ Gubser and Peebles [15] 
considered structure formation in the presence of both scalar and gauge long-range forces, but concentrated on a region 
of parameter space in which the gauge fields were subdominant. Refs. [TOJIT7] considered several models for a hidden 
dark sector, including one manifestation in which the dark matter consists of heavy hidden-sector staus interacting 
via a copy of electromagnetism. The effect of dimcnsion-6 operators containing a new U(l) gauge boson and SM 
fields was considered in Ref. |18j . for models where the only fields in a hidden sector are charged under the new force. 
Additional models which contain unbroken abelian gauge groups may be found, for example in Refs. |19l 120] . In this 
paper, we construct minimal models of dark matter coupled to a new unbroken U(l) £>, leaving the dark fine-structure 
constant and dark-matter mass as free parameters, and explore what regions of parameter space are consistent with 
astrophysical observations and what new phenomena might arise via the long-range gauge interaction. 

In Section [IT] we introduce our notation for a minimal dark-matter sector including a new abelian symmetry U(1)d- 
We then consider the bounds on the new dark parameters from successful thermal production of sufficient quantities of 
dark matter as well as requiring that BBN and cosmic microwave background (CMB) predictions remain unchanged. 
The restrictions of parameter space are closely related to those resulting from standard short-range WIMP dark 
matter. In Section |III[ we consider the effect of long range interactions on DM particle interactions in the halos 
of galaxies. By requiring that our model not deviate too greatly from the predictions of collisionless DM, we find 
that the allowed regions of a/m x parameter space from Section UU are essentially excluded. In order to evade these 
constraints, Section IV describes an extended model, where the dark-matter candidate is charged under both SU (2) j, 
and the new U(1)d- Additional effects of dark radiation are presented in Section |VJ and we conclude in Section VI 



We note that our model does not address the hierarchy problem, nor provide a high-energy completion to the SM. 
However, new gauge groups and hidden sectors may be generic results of many such high-energy theories (e.g. string 
and grand unified theories), and a WIMP coupled to an unbroken U(l) is certainly a plausible low-energy manifestation 
of such theories. The most important lesson of our model is that interesting physics might be lurking in the dark 
sector, and it is worthwhile to consider a variety of possible models and explore their consequences for astrophysics 
and particle physics. 

II. DARK RADIATION AND THE EARLY UNIVERSE 

We postulate a new "dark" abelian gauge group U(1)d with gauge coupling constant g and dark fine-structure 
constant a = g 2 /An. In the simplest case, the dark matter sector consists of a single particle x with U(l) d charge of +1 
and mass m x along with its antiparticle \. For definiteness, we take \ to be a fermion, though our results are applicable 
to scalars as well. As the limits on new long range forces on SM fields are very stringent, we assume that all the SM 
fields are neutral under U(1)d- For the moment we take the x ne ld to be a singlet under SU(3)c x SU (2) l x U(1)y, a 
restriction that will be relaxed in Section [rV] As a result, this extension of the SM is anomaly free. In this Section, we 
will derive constraints on the mass m x and coupling a from the evolution of dark matter in the early universe. Two 
considerations drive these constraints: the dark matter must provide the right relic abundance at thermal freeze-out, 
and the dark radiation from the U(1)d cannot contribute too greatly to relativistic degrees of freedom at BBN (a 
similar bound coming from the CMB also applies but is weaker). 

The degrees of freedom in the dark sector are thus the heavy DM fermions x an d massless dark photons 7. The 
Lagrangian for the dark sector is 

C = x(0 + m x ) X -\p flv F^. (1) 

Here = 9 M — igA^ and F^ v is the field-strength tensor for the dark photons. We assume that the mixing term 
cF^ iy F^ v is set to zero at some high scale (say the GUT scale). This is a self-consistent choice, since if there is no 



mixing between the dark and visible sectors, c = is preserved by the renormalization group evolution. (In Section IV 
we argue that mixing is not generated by radiative corrections even when x carries SU(2)l quantum numbers.) This 
choice allows us to bypass constraints on a new U(l) coming from mixing between the photon and dark photon, that 
is, 'paraphotons' 122j . We have no a priori assumptions on the parameters m x and a, though as we shall see, it 
suffices to think of the former as 0(100 - 1000 GeV) and the latter < O(10" 2 ). 

We now follow the thermal history of the dark sector. Our analysis follows that of Ref. [17]; we rehearse it in a 
slightly simpler context here to illustrate how the results depend on our various assumptions. If the visible sector and 



1 Broken U(l) forces have, of course, also been considered, see e.g. Ref. [9] 
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the dark sector are decoupled from each other, they may have different temperatures, T and T, respectively; rapid 
interactions between them would equilibrate these two values. After inflation, the two sectors could conceivably reheat 
to different temperatures, depending on the coupling of the inflaton to the various fields. Even if the temperatures 
are initially equal, once the two sectors decouple as the universe expands and cools, entropy deposited from frozen- 
out degrees of freedom in one sector will generally prevent the dark temperature T from tracking the visible sector 
temperature T. The ratio 

f = f/T (2) 

will depend on the spectrum of both sectors, and is itself a function of T. 

As the temperature drops below a particle's mass, the associated degrees of freedom freeze out and dump entropy 
into their respective sectors (dark or visible). This causes the temperature of that sector to decline more slowly than 
1/a, where a is the scale factor of the universe. As the entropy density s of the visible sector and s of the dark sector 
are individually conserved after decoupling, we must separately count the degrees of freedom in these two sectors. 
There are two definitions of degrees of freedom of interest to us: g„ and g*s- The former is defined as 



E 



t) + 8 E 

i— fcrmions 



(3) 



and is used in calculation of the total relativistic energy density, pji oc <?*T 4 . Here <?, is the number of degree of 
freedom for particle species i, Ti is the temperature of the thermal bath of species i, and T is the temperature of the 
photon bath. The sums run over all active degrees of freedom at temperature T. Separating out the visible fields, 
can be written as 



E + 1 E (4) 



5*vis ~f 

i— bosons i— fcrmions 



where the sums now run over the dark particles. If we restrict the visible sector to the SM, then the term <7* v is is 
106.75 above the top mass, dropping gradually to ~ 60 at T = Aqcd- Between 100 MeV > T > 1 MeV, g* v is = 10.75, 
and drops again to 3.36 in the present day. (See e.g. Ref. [21] for more detail.) 

Similarly, the total entropy density s to t ( a conserved quantity) at a photon temperature T is proportional to g^T 3 , 
where 



3 - /m \ 3 



E +1 E MS) (*) 



K T J 8 ^ ° \T 

i— bosons i— fcrmions 

= 3*s,vis + £ g^{Tf + 7 - g£{Tf . (6) 

i— bosons i— fcrmions 

Prior to neutrino decoupling, all the relativistic standard model degrees of freedom are in thermal equilibrium at a 
common temperature. Thus, before T ~ 1 MeV when neutrinos decouple, W6 have £/*vis — *?*S,vis- 

Furthermore, we 

may split the dark g*§ into heavy and light degrees of freedom: (/heavy and (/light, where the heavy degrees of freedom 
are non-relativistic at BBN. We are interested in the number of degrees of freedom at BBN (T ~ 1 MeV) because 
formation of the experimentally observed ratios of nuclei are very sensitive to the expansion of the universe at that 
time, which is related to the energy density of radiation through the Friedmann equation. From this, a bound on the 
number of relativistic degrees of freedom can be derived |17j . 

Using the separate conservation of the visible and dark sector entropy and the previous definitions, we see that, at 
BBN 



Slight£(7BBN) _ 9*vis(TbBn) 



(ffheavy + Slight ) £(TrH ) 3 <7*vis (TrH ) 



(7) 



here we have set g*s,vis = 5*vis (recall that g* v is(2BBN) = 10.75). The BBN bound on relativistic degrees of freedom 
is usually stated in terms of number of light neutrino species in thermal equilibrium at the time: N u — 3.24 ±1.2 
[2"5] . Here the error bars correspond to 2a (95% confidence). Therefore, assuming three light neutrino species in the 
visible sector, if the dark sector is not to violate this bound, we must require 

5ii g htC(?BBN) 4 = -x2x(JV„-3)< 2.52 (95% confidence). (8) 
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Combining Eqs. ^ and we find that 



Slight 



ffheavy + flight 10-75 
Slight S*vis (2hh) 



4/3 



£(7rh) 4 < 2.52 (95% confidence). 



(9) 



Since the high energy completion of the visible sector must at minimum include the SM fields, SWsC^rh) > 106.75; 
a bound on the dark sector flight and <?hcavy can be derived for a fixed value of £(7rh) (see Fig. lj). Increasing the 
number of visible sector degrees of freedom at high temperatures (for example to that of the MSSM) relaxes this 
bound. 





FIG. 1: The allowed values of dark <?n g ht (those degrees of freedom relativistic at Tbbn ) and <?hcavy (the remaining dark 
degrees of freedom) arising from BBN constraints Eqs. |8| and |9j. The allowed regions correspond to 95% confidence levels for 
£(7rh) = 1 and a visible sector g tv i s — 106.75 (red), f (Trh) = 1 and g* v i B — 228.75 (corresponding to MSSM particle content, 
in blue), and £(Trh) = 1.4(1.7) and gr„ v is = 106.75(228.75) (in yellow). The minimal dark sector model of this paper is noted 
by a black star at <?n g ht = 2 and <?hcavy = 3.5. 



In the case of £(Trh) = 1, we see that the minimal model of the dark sector (only heavy xlx an( i light 7) is safely 
included. Due to the fourth power of £ entering into Eq. (J9|, if the minimal dark sector is not to be ruled out, we 
find £(Trh) < 1.4(1.7) for the SM(MSSM) particle content. A similar bound on relativistic degrees of freedom can be 
derived from the cosmic microwave background, but provides a weaker la exclusion limit |171 124j . 

We now turn to bounds on the coupling a and dark matter mass m x coming from the dark matter abundance. 

At temperatures T much above m x , the x particles are kept in thermal equilibrium with the dark photons 7 (and 
possibly other particles in the dark sector) via pair annihilation/creation as in the Feynman diagrams of Fig. [2j Since 
the annihilation can proceed via s-wave processes, the thermally averaged cross section (av) is, to leading order, 
independent of v: 2 

(av)*a = ^+O(v 2 ). (10) 

Using this, the relic density of the x particles may be easily calculated (see, for example Ref. [21]). 

As a rule of thumb, the dark matter drops out of thermal equilibrium when the rate T of annihilation XX ~> 77 (and 
the reverse process) is outpaced by the expansion of the universe H . Using the Boltzmann equation, the contribution 
of x to the energy density of the universe can be more precisely calculated as 

An+l)x n f +1 GeV- 1 
n D Mfr 2 = 1.07xl0 9 — ^— . (11) 



2 Strictly speaking, there will be a Sommerfeld enhancement in this cross section in the limit v — > 1251 . This will slightly change the 
relic abundance |26j . but we leave the detailed analysis for future work. 
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FIG. 2: Pair annihilation/creation of dark matter \ m t° dark photons 7 via f and 11-channel exchange diagrams. These 
processes keep the dark sector in thermal equilibrium until the % particles become non-relativistic. 



Here Xf is the ratio m x /Tf where Tf is the dark temperature at time of freeze-out and n — for s-wave processes. 
The quantity Xf is given by 



In 



0.038(n + l) 



mp\m x <7o 



In In 



0.038(n+ 1) 



'9* 



(12) 



where g is the number of degrees of freedom in the \ system (namely 4) . 

As enters into the formula for Xf only logarithmically, we may make the approximation that g*s ~ 100 if x 
freezes out while T is above Aqcd ■ We make the additional assumptions that the only degrees of freedom in addition 
to the SM are the 7 and \ m the dark sector and that £(Trh) = 1. We shall consider how these assumptions may be 
relaxed later. 

Under these assumptions, the contribution of the dark sector to g* and g*s is 2+ (7/8) x4 = 11/2. As no dark degrees 

1 /3 

of freedom have frozen out yet, £(T f ) = ( g 9 "' f^) ) ?( t Rh) ~ 1- With the measured value tt DM h 2 = 0.106 ± 0.08 

[1], we may solve for the allowed values of a as a function of m x in Eq. (111. The resulting band is shown in Fig. [3j 

In this discussion we have assumed that the process which sets the relic abundance of x is annihilation into 7s, as 
shown in Figure [2] As we will argue in the next section (and as is already shown in Figure j3|, the values we obtain 
for a from this calculation are incompatible with bounds from galactic dynamics unless m x > 10 5 GeV (at which 
point a is non-perturbative). However, we can get the correct relic abundance even with much lower values of a by 
adding other annihilation channels, such as the weak interactions, as explored in Section |IV| In that case, the "relic 
abundance allowed region" discussed here really becomes an upper limit; if the dark fine-structure constant is larger 
than that value, annihilations are too efficient, and the correct abundance cannot be obtained. 

We now consider how changing our assumptions on g* and £ can change our conclusions on the allowed parameter 
space. The parameter £(T/) does not enter explicitly into the calculation for JIdm h 2 , however it does affect the 
number of active degrees of freedom at freeze-out directly, through Eqs. Q and ([6|, and indirectly by allowing the 
temperature T to differ from T. If £ < 1 , T < T and there could be many more heavy visible degrees of freedom still 
active when x freezes out. £ > 1 would reduce the visible degrees of freedom. However, as we have seen in Eq. (|9j, it 
is difficult to construct a scenario with large £, short of a massive increase in g* V is and small values of <?* heavy +5*iight- 
We include in Fig. [3] the bounds from both a large and small value of <7*. The large limit is g* v is(Tf) = 228.75, 
(i.e. equivalent to the MSSM degrees of freedom), £(Trh) = 1, and (/heavy + Slight = 100, while the small value is given 
by 9*vis{Tf) — 60, (i.e. equivalent to the SM degrees of freedom at Aqcd)i £(?rh) — 0.1, and $hcavy + <?iight — 5.5. 



III. GALACTIC DYNAMICS 



Although freezeout in our scenario is similar to that in the standard WIMP scenario, the long-range DM-DM 
interactions implied by the unbroken U(1)d may lead to considerably different DM phenomenology in the current 
Universe, and in particular in galactic halos. In this scenario, dark-matter halos are composed of an equal mixture of 
X and x- The overall halo will be U(1)d neutral, eliminating long-range forces that are incompatible with experiment. 

However, nearest-neighbor interactions between x particles remain, and these interactions can be constrained by 
observations that suggest that dark matter is effectively collisionless. Constraints to dark-matter self-interactions arise 
from evidence for nonsphcrical cores for some dark-matter halos (collisions tend to make the cores of halos round) 
[27] and from evidence for dark-matter halos with large phase-space densities (collisions would reduce phase-space 
densities) 7, 28, 29 . Roughly speaking, a bound to DM-DM interactions can be derived by demanding that scattering 
induces no more than a small fractional change in the energy of a typical DM particle in a galactic halo during the 
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m x (GeV) 

FIG. 3: The allowed regions of d vs. m x parameter space. The relic abundance allowed region applies to models in which 
U(1)d is the only force coupled to the dark matter; in models where the DM is also weakly interacting, this provides only an 
upper limit on a. The thin yellow line is the allowed region from correct relic abundance assuming Odm/i 2 = 0.106 ±0.08, 
£(7rh) = 1, <?*vis ~ 100, and ghcavy + ffiight = 5.5 while the surrounding blue region is g tv i s — 228.75(60), £(Trh) = 1(0.1), and 
ghoavy + Slight = 100(5.5) at the lower(upper) edge. The diagonal green line is the upper limit on a from effects of hard scattering 
on galactic dynamics; in the red region, even soft scatterings do not appreciably affect the DM dynamics. We consider this to 
be the allowed region of parameter space. 



history of the Universe [BJ. This translates to an upper bound of ~ 0.1 cm 2 /g on the more familiar quantity a/m x . 3 
A separate bound of cr/m x < 1.25 can be derived from the Bullet Cluster [301 131], but as this is less restrictive we 
ignore it here. 

To illustrate, we first consider hard scattering of a x off another \ or Xi where energy on the order of m x v 2 /2 is 
exchanged. The mean free time t for a x to undergo a hard scattering with another x(x) is given by 

T=T- i —y (13) 

where n is the number density of dark matter, a is the hard-scattering cross section, and v is the velocity of the 
dark-matter particles. The number N of dark-matter particles in the Galaxy is 



N = Mcal « 10 64 
m x 



and n ~ 3iV/47ri? 3 , where R is the radius of the Galaxy. The velocity v is 



GM Gal GNm x 

^^"V^^- (15) 

The dynamical time Td yn in the Galaxy is 

r dyn = 2irR/v. (16) 



3 This can be seen from Eq. (131 g using the age of the universe for t, and Galactic parameters p = nm x = 0.3 GeV/cm 3 , v/c = 10 3 . 
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Taking Td yn ~ 2 x 10 8 years for the Milky Way, the average time for a hard scatter for a dark-matter particle is greater 
than the age of the universe if 

r dyn 3Na ~ y ' 

A hard scatter occurs when two particles pass close enough so that their kinetic energy is comparable to their 
potential energy. The impact parameter that defines a hard scatter is thus 

2d 

iWd=-~ ■ (18) 

v z m x 



Taking the cross section for hard scatters to be Uhard ~ ^hardi an d using Eq. (15 1 for v, we find 



Twd G 2 miN ^ 

— = * > 50. (19) 
r dyn 6a' 

Using G = nip 2 w 10 -32 TeV -2 we find the hard scattering limit on the U(1)d coupling constant to be 

The allowed region arising from this bound is shown in Fig. [3] 

We now turn to the effect of soft scattering on the allowed values of a and m x . Here we consider the approach of 
one x particle towards another x(x) at impact parameter b. By definition, for soft scattering b > fehard- The velocity 
change induced by the encounter is 

5v = ±^* r . (21) 
m x bv 

As one dark-matter particle orbits the Galaxy, it sees a surface density N/irR 2 of dark matter. The number of 
interactions that occur between an impact parameter b and db is Sn = (N/TrR 2 )2irbdb. While the change in 5v over 
these interactions should average to zero, this is not true for 6v 2 : 

Sv 2 = (dvfdn = ^fL b-i-db. (22) 

Integrating Sv 2 from &hard to the maximum possible impact parameter in the Galaxy, R, gives the total change in v 2 
as the particle orbits once through the halo: 



ia 2 N 1 . „ . 8a 2 N , ( GNm 2 x 



AV = 2 2m MR/Kzrd) = a 2 R 2 ln ' ( 23 ) 

miv z R z miv z R z \ 2a j 



x x 



The number r/rd yn of orbits it will take for the dark-matter particle to have Av 2 /v 2 ~ 0(1) is 

(24) 



Tsoft _ G 2 m\N ]n _ 1 (GNm 2 x \ ^ 

Td yn 8d 2 I 2fl J ~ 



The logarithmic suppression in Eq. ( |24| relative to Eq. (19 1 is due to the long-range Coulomb force generated by 
the U(1)d- As can be seen in Fig. [3] the allowed region from these considerations of Galactic dynamics completely 
exclude the a/m x band that gives the correct relic abundance up to m x ~ 30 TeV. For m x ~ 1 TeV a dark matter 
candidate which freezes out due to U(V)d interactions is ruled out from such considerations. In particular, models 
such as that in Ref. |17j with m x ~ mw and a hidden copy of electromagnetism (i.e. a — 1/137) are ruled out, even 
though the freeze-out proceeds through hidden-sector weak interactions rather than a U(1)e>. Interestingly, a = a is 
allowed for m x > 2 TeV. 

Before considering whether such a model may be valid if our assumptions are loosened, we should ask why Galac- 
tic dynamics do not similarly exclude WIMP dark matter. After all, both models have similar cross sections for 
annihilations in the early universe (Eq. ( |10[ )) as is required for the correct relic density. Though the soft scattering 
limit clearly will not apply due to the short range nature of the broken SU(2) l, naively it would seem that the hard 
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scattering limit Eq. (|19j) should apply to WIMPs equally well. However, notice that the threshold for hard scattering 
with a U(1)d is dependent on energy. As the temperature drops, the cross section rises, as the x particles no longer 
have to approach as close in order for to potential energy V(r) to be of the order of the kinetic energy. Contrast this 
to hard scattering from WIMPs, where the cross section is always proportional to a 2 /m'^ )M , regardless of the velocity. 



Entering this cross section into Eq. (17 1, results in the uninteresting bound that ttidm ^5 10 13 TeV for WIMP dark 
matter from Galactic dynamics constraints. 

It is difficult to see any way of avoiding the bounds from Galactic dynamics, so we look to loosen the limits derived 
in Section |nj Clearly if the interaction responsible for freezing out the relic density is not the U(1)d constrained by 
soft scattering, then d < 10~ 3 is not ruled out. We consider such examples in the next Section. However, we first 
consider the possibility that our assumptions in deriving the relic density are too conservative. 

From Eq. ( fTT| , if we reduce d (and therefore a ) in order to satisfy the scattering bounds, we must either decrease 
x f or increase g*s/\/9*- In lowering d by a minimum of two orders of magnitude, x / /(#*s I yfs*) must likewise increase. 
As Xf depends only logarithmically on d and the number of degrees of freedom, so it is hard to see how it it could 
be increased sufficiently to counterbalance d of order 10~ 3 (rather than d ~ 1CP 2 ). We conclude that the number 
of effective degrees of freedom must be increased. From Eqs. Q and (|6j, we see that if £ = 1, then at freeze-out we 
must have 



, E + 5 E .'/<-- |()J - (25) 



2=fermions 

From Eq. these ~ 10 4 degrees of freedom must exist in the visible sector at Tf, rather than the dark sector. 

Alternatively, we could imagine that there are no (or few) new particles beyond the minimum x an d 7 at freeze-out, 
but instead £ 3> 1 . In this limit 

^«£~ 10 2 . (26) 

This limit is more troublesome; from Eq. (|J we saw that large values of £ at the reheating scale (and subsequently 
Tf) very quickly violate the bounds on relativistic degrees of freedom at BBN. Clearly, by increasing the number of 
degrees of freedom in the visible sector, this bound could be avoided as well. However, we are left with the conclusion 
that either £(7/) ~ 10 2 or there exist ~ 10 4 new particles at a few hundred GeV to a TeV. We leave it to the reader 
to decide how palatable these alternatives are. 

A separate, but conceptually similar, bound on scattering can be placed by considering the interaction of galactic 
dark matter with the hotter DM of the surrounding cluster. Scattering will cause heating in galactic DM, and 
eventually evaporate the halo. From Ref. [32] the characteristic time for this evaporation is given by 

W = 3.5 x 10 9 years ( ^fV {^^rY ff?" X' ■ (27) 

\cm /g/ \lCrkm/s/ \ 1.3 x lO^Mopc - - 3 / 

We may estimate the cross-section for soft-scattering by calling the path length A over which a single particle loses 



of order its initial kinetic energy (Av 2 /v 2 )~ 1 R, where R is the radius of the galaxy, and Av 2 /v 2 from Eq. (23 1 is the 
fractional energy loss as the particle travels once through the halo. This can be expressed as an effective scattering 
cross section by setting A = (ner) -1 , where n = N/R 3 is the number density of DM in the halo, we find 

R™ 2 / GNm 2 \ 

In ■ (28) 



m x m 3 i> 4 I 2d 

Letting the cluster velocity and density take on the canonical values (^cluster — 10 3 km/s and p c iustcr — 1-3 x 
1CP 3 Mqpc -3 , where Mq is the solar mass), we can place limits on d and m x by requiring that i cva p. is longer 
than the age of the universe. Numerically, we find this bound less stringent than that from soft-scattering of particles 
within the Galactic halo, Eq. |24} 

It is interesting to note that, aside from logarithmic enhancements, the bound placed on d vs. m x parameter space 
from soft scattering is essentially a line of constant cr/m x (that is, they are, up to log corrections, lines of slope 2/3 
on the log-log plot). As mentioned, limiting DM to one hard scattering in the lifetime of the universe is equivalent to 
bounding cr/m x in the Galaxy to be ^0.1 cm 2 /g. It has been suggested in the literature that values of c/m x in the 
range 0.01 — 5 cm 2 /g [71 H?l I2H1 US] may provide better agreement between simulation and observation. Therefore, our 
limit from soft-scatterin g sh ould be considered as the general region at which interaction effects may become relevant. 



Additionally, from Eq. (28 1 as cr/m x cx v~ , it should be expected that the soft-scattering bound will vary greatly 
in DM systems with a range of virial velocities v. In particular, we surmise that a bound even stronger than that 
estimated here can be obtained from the dwarf galaxies that exhibit the highest observed dark-matter phase-space 
densities 133 . 
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IV. WEAKLY COUPLED MODELS 



In this Section, we examine an expanded version of our minimal model: one in which the x dark matter particles 
possess SU (2) l quantum numbers in addition to a U( 1) d charge. For such SU (2) i x U(l) d particles, the cross section 
for freeze-out in the early universe is dominated by the weak interaction a ~ a 2 /m x , and the U(1)d contribution 
is negligible for the small values of a under consideration. At late times the situation is reversed. The weak cross 
section remains small, as it is the result of a short-range force. Howe ver the long range cross section for soft scattering 
increases as the dark matter cools and slows, as exemplified in Eq. (23 1. This allows the strength of a to be ~ 1CP 3 
as required by Galactic dynamics without running afoul of the relic density conditions, which would require a ~ 1CP 2 
(when m x ~ 1 TeV). 

We therefore take our Dirac fermion \ to be a (1, vl)y,d multiplet of SU (3)c x SU(2) l X U(l)y x U(1)d, where we 
shall take the U(1)d coupling to be in the region of Fig. [3] allowed by soft scattering. Thus a < 1CP 3 . The behavior 
of this model in the early universe is very similar to the 'minimal dark model' of Ref. [34 , from which we take many 
of our constraints. 

In outlining our original model in Section [n] we set the coefficient of the mixing term F^ V F^ U to zero at the high 
scale. Clearly loops involving \ would generate a non-zero mixing if the \ field possesses non-zero hypercharge Y . In 
order to avoid this complication, we set Y = 0. 

Our x particle must be neutral under U(\)em- With the assumption of Y = 0, this requires \ to sit in an n-plet 
of SU{2)l where n is odd (i.e. n = 3, 5, . . .). In the spirit of simplicity we take n = 3, so the \ triplet contains the 
neutral x° an d (electromagnetically) charged x^i a U with U(1)d charges of +1. Due to SU(2)l loops, the x are 
166 MeV heavier than the x°i an d decay before BBN. If the dark matter mass is m x = 2.4 TeV, the correct dark 
matter abundance (including production and then decay of x ± ) results from thermal freeze out (see Ref. [34 )■ We 
note that our model does have the nice feature of automatically suppressing unwanted decays of x m to SM particles, 
as by assumption x is the lightest particle charged under U(V)e>. 

This minimal model is anomaly free. Triangle diagrams with one or three SU (2) l vertexes vanish by the tracelcssncss 
of the SU (2) l generators. The diagrams consisting of an odd number of U(l) d vertexes also vanish as the dark sector 
contains only two Weyl fermions, one with +1 under U(l)r>, and the other with —1. 

This model does not run afoul of BBN (or CMB) bounds. As in the pure U(1)d theory, the only new relativistic 
degrees of freedom at BBN are the two from the 7. Due to the interactions between x an d the weakly charged SM 
fields, we expect the temperatures T and T to track, so £ = 1 until the x freeze-out. With small values of d, the dark 
photons may freeze-out earlier, and would thus be colder. However, if we take the worse-case scenario that the dark 
photons do not decouple until after the x undergo freeze-out we find (from Eq. ^ ) that BBN bounds are satisfied as 
long as freeze-out occurs when 



9* 



> 18.8. 



(29) 



This is easily satisfied for any model that freezes out before the QCD phase transition. 

Next we must check that our x does not have too large of a coupling to SM particles. We first demonstrate that 
no mixing occurs between the photon and the dark photon 7. As indicated previously, we assume that there is no 
F iLV F lil ' term at high energies. With purely SU(2)l x U(1)d coupling, we find that the diagram Fig. [4^, vanishes. 
This is because any such vertex can be rewritten as the 7 coupling to a % or % which then couples to the 7 through 
some vertex involving SM fermions and SU(2)l couplings (Fig. |4|d). However, since the mass and SU(2)l couplings 
of x are the same as those of x Y e t the U(1)d charge is opposite, the sum of the two diagrams is zero. 





FIG. 4: Feynman diagrams leading to 7/7 mixing. The vertex in a) can be expanded into that shown in b), as the only particle 
to which the 7 couples is x/X- Since the mass and SU(2)l charge of these two particles are the same, yet they possess opposite 
U(1)d charge, the sum of the \ and % diagrams in b) is zero, and the overall mixing vanishes. 



Similarly, the coupling between 7 and a standard model fermion / is also zero. The relevant diagrams are shown 
in Fig. [5] Again, the vertex between / and 7 (Fig. [5^,) can be divided into the x/x vertex connecting with 7 and a 
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vertex between x/x vertex connecting with / (Fig. |5p). As the latter vertex is identical for x an d X but the former 
has opposite signs, the overall diagram vanishes. 



a) 

/ 




b) 

/ 


xlx 














/ 




/ 







FIG. 5: Feynman diagram leading to 7 interactions with SM fermions /. The vertex in a) can be expanded into that shown in 
b), as the only particle with an interaction with 7 is the xlx.- Since the mass and SU(2)l charge of these two particles are the 
same, yet the U(1)d charges are opposite, the sum of the x an d X diagrams in b) is zero, and the overall coupling of / to 7 is 
therefore zero as well. 

The lowest order coupling of SM fermions to 7 occurs at a 2 a. This is due to a two loop effect, as shown in Fig. [6] 
and unlikely to be accessible in direct detection. We can represent this interaction by an effective Lagrangian whose 
lowest order term is given by t^-F^F^" ff where p = Xf 9 ^ and A/ is the Yukawa coupling of the fcrmion that is 

involved. Let us estimate the order of magnitude of this interaction. To be conservative we use the Yukawa coupling 
of a u quark and take a = 10~ 2 ; which by galactic dynamics is the maximum allowed value for m x ~ 2 TeV. With 
these values we find j3 ~ 10~ 10 and J^- ~ 10~ 20 GeV~ 3 . We estimate that the interaction length for dark photons 

inside the cores of stars would be on the order of 10 18 km, and thus this interaction would not introduce a potentially 
dangerous new source of stellar cooling. 









7 




7 




xlx 




w a <> 






f f 



FIG. 6: The leading order interaction of the dark sector with SM fermions. The dark photons 7 couple to a loop of x particles, 
which couple through two SU(2)l gauge bosons to SM fermions. Coupling through a single SU(2)l boson is zero due to the 
tracelessness of r a . 

Due to the high-order interaction between 7 and SM particles, we cannot expect to directly observe the dark 
radiation. In addition, while the x fields would have a direct detection cross section of 10 -44 — 10~ 45 cm 2 |Mj and 
so could be seen in SuperCDMS, any such detection would be indistinguishable from a scenario without the dark 
photons. Therefore, the presence of a new unbroken U(X)u in the dark sector could only be probed via its effect on 
galactic dynamics. Clearly in the limit that a — > 0, the Galactic structure would remain unchanged. Values of a near 
the maximum allowed from soft-scattering (i.e. a ~ 10~ 2 for the SU(2)l triplet candidate with m x ~ 2 TeV) should 
have a measurable effect on the halo structure, as in this regime the dark matter is no longer completely collisionless. 
A full study of this effect requires simulations beyond the scope of this paper, though some additional considerations 
are discussed in the following Section. 
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V. OTHER EFFECTS OF DARK PHOTONS 

The existence of a dark matter 'plasma' may have additional effects that could significantly affect structure for- 
mation. We mention three possibilities here: bremsstrahlung, early universe structure formation, and the Weibel 
instability in galactic halos. The first two result in much weaker bounds than those already derived, and are men- 
tioned here only for completeness. The Weibel instability may have significant and visible effects in the halo, but 
requires simulation beyond the scope of this paper. 

A. Bremsstrahlung 

Emission of a soft 7 during a x/x collision could conceivably serve as another energy loss mechanism in the halo on 



par with soft and hard scattering as outlined in Section III To derive a bound on a as related to m x , we make the 
same assumption as in the case of soft scatter: over the lifetime of the universe, a dark-matter particle cannot lose 
on order of its initial kinetic energy through bremsstrahlung of dark radiation. By assuming dipole radiation during 
a soft collision, we find that 



3 Gm 3 x R _ x ( GNrn 2 ^ 
64 a 3 n I 2a 



> 50. (30) 
However this bound is weaker than that from both hard and soft scattering over the parameter space of interest. 



B. Structure Formation 

In the early universe, structure cannot grow until after matter/radiation equality. Until the matter (which can 
clump) decouples from the dark radiation (which cannot), density perturbations remain fixed. We can estimate the 
scale factor at which this occurs by finding the redshift z + at which the dissipation time (the time over which the 
velocity of a dark matter particle is significantly perturbed by the radiation) becomes longer than the Hubble time 
if -1 . The argument follows that in Ref. [55] for the decoupling of baryons from the photon bath. 

The dissipation time is the logarithmic derivative of the velocity: 



-1 _ „.-i^ _ ..-1 F 
dt 



Mb.* = ' "77 - ' ~~ ' ( 31 ) 



Here F is the force due to radiation pressure, 

4 



F=^a T aT 4 v, (32) 



where 

8tt a 2 
3 m 



°t=-^ (33) 



x 

is the Thomson cross-section for dark matter interacting with dark photons and (as before) T is the temperature of 
the dark photons. As we shall see, the decoupling occurs when the universe is radiation dominated, so the Hubble 
time is given by 

4tt 3 T 4 

" 2 = < 34 > 

Here T is the photon temperature. 

The conservation of entropy relates the photon temperature T at redshift with the photon temperature today, 
To, 

T=( 9 -^4) 1/3j °. (35) 
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Combining Eqs. (31) and (34 1, we find the decoupling redshift to be 



1 + z* 



3 

Tg 



& 2 Tq mpi 



9*{T) 



1/2 



18 e-4 



2.3 x 10 1S £ 



10" 



n 



TeV 



9*s(T ) 



2/3 



1/2 



9*s(T) 
9*s(T ) 



2/3 



(36) 



As before £ is the ratio of dark photon temperature to photon temperature at redshift (recall that it is difficult 
to construct models where £ is much larger than unity) . The number of degrees of freedom that contribute to the 
entropy density today, 5 (To), is of order unity. The decoupling occurs extremely early, before even dark matter 
freeze-out. 4 As a result, it seems that this effect will be cosmologically irrelevant. 



C. Plasma Instabilities 



In Section III we constrained a by demanding that dark matter be effectively collisionless in galactic halos, under 
two-body interactions. However, there may be collective plasma effects that affect DM dynamics on timescales much 
shorter than those due to two-body interactions. Unfortunately, it is difficult to state with confidence what the 
observational consequences of those effects will actually be, even if they are relevant. Given theoretical uncertainties 
about the nonlinear evolution of such instabilities, we leave the detailed implications to future work. 

As a simple example we consider the Weibel instability |36j . an exponential magnetic-field amplification that arises 
if the plasma particles have an anisotropic velocity distribution. Such anisotropies could arise, for example, during 
hierarchical structure formation as subhalos merge to form more massive halos. Similar instabilities in the baryonic 
gas have been postulated to account for the magnetic fields in galaxy clusters [37] . The growth rate T of the magnetic 
field is 



c 



where uj p is the plasma frequency, p sa 0.4 GeV/cm 3 is the dark-matter density, and v is the velocity of the dark 

1CT 3 , > 

r-itrV 1 x " (38) 

(m x /TeV) 

To be relevant for galactic-halo formation, the timescalc for magnetic- field amplification should be shorter than 
the dynamical timescale t of the merging subhalos. The instability will be therefore be of interest when 

(&) 5 10,141/5 fe) • (39 > 

This range of a and m x encompasses the entire parameter space of interest for any reasonable value of r. Therefore, 
we suspect that galactic structure will be affected by plasma effects in the dark matter due to the U(1)d even when 
a is not near the boundary of allowed values from soft scattering. One possibility is that nonlinear evolution would 
result in a strongly magnetized plasma, and if so, dark matter would be effectively collisional and thus probably 
inconsistent with data. However, theory and simulations that study the nonlinear evolution of the Weibel instability 
for relativistic pair plasmas and nonrelativistic electron-proton plasmas do not yet agree whether the magnetic fields 
survive, and simulations for the equal-mass nonrelativistic plasma we are considering have not been performed. It 
is therefore premature to conclude that these instabilities will result in effectively collisional dark matter; a more 
detailed study will be required to assess these effects. 




= \ o ' ( 37 ) 



matter within the colliding halos. Assuming v/c ~ 10 3 , we find 

a 1 ' 2 



4 This is not a contradiction: freeze-out is the time when the dark particles and antiparticles stop annihilating, while decoupling occurs 
when the dark photons stop imparting significant velocity to the dark matter. 
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VI. CONCLUSIONS 

Given how little direct information we have about the nature of dark matter, it is of crucial importance to explore 
models in which the DM sector has an interesting phenomenology of its own. In many ways, an unbroken £7(1) gauge 
field coupled to dark matter is a natural way to obtain a long-range interaction between DM particles. In contrast to 
the case of hypothetical long-range scalar fields, the masslessness of the gauge field is protected by a symmetry, and 
the absence of long-range violations of the equivalence principle is naturally explained by the overall charge neutrality 
of the dark plasma. New unbroken £7(1) 's can appear naturally in unified models. 

While a dark U(l) may be realized as a broken symmetry with massive vector bosons, it has been pointed out that 
there are few constraints on the massless, unbroken case from the early universe. We have verified that the minimal 
model, with just a single massive Dirac fermion for the dark matter and a massless dark photon, is consistent with 
limits obtained from the number of relativistic degrees of freedom at BBN, with relatively mild assumptions on the 
reheating temperature of the dark sector. More complicated models are also allowed, depending on the details of 
spectrum and reheating. 

We found that one cannot build a dark matter model charged under a hidden unbroken U(1)d m which this new 
gauge group is responsible for thermal freeze out. As can be seen in Fig. [3j the required values of a and m x required 
for the x particles to form a thermal relic would violate bounds coming from limits on hard and soft scattering of dark 
matter in the Galactic halo. As an important consequence of this argument, models in which dark matter couples to 
an exact copy of ordinary electromagnetism (in particular, with a = a) are ruled out unless m x > a few TeV. This 
constrains the parameter space of models with hidden copies of the SM or the MSSM in which the dark matter is 
electrically charged, such as the model in Ref. [T7] where the stau was suggested as a dark matter candidate. 

By adding additional interactions to increase the annihilation cross-section, it is possible to build a scenario with 
an unbroken dark U(l) and the correct relic abundance. Introducing another short-range force coupling to the \i f° r 
example the familiar SU{2)l, can provide an appropriately large cross section for x/x annihilation. The new coupling 
a must then be relatively small (compared to the SU(2)l oc) in order to evade Galactic dynamics bounds. 

The simplest model which realizes this situation is a Dirac fermion in a triplet of SU(2)l (in order to avoid 
U(1)y /U(1)d mixing). Bounds from the early universe then force m x to be on the order of a few TeV, which implies 
a < 10 -2 . Since all couplings between the dark radiation and the SM enter at two loops (and require two dark photons 
in the process), it would be very difficult to observe the presence of the new gauge group through direct detection. 
Instead, the best search strategy would be an indirect one: looking for the effects on Galactic dynamics arising from 
a soft scattering mediated by a long-range force. Clearly, as a goes to zero, the model becomes indistinguishable 
from minimal weakly coupled dark matter. However, if the coupling is near the limit from soft scattering, one would 
expect detectable deviations from the assumptions of collisionless dark matter currently used in simulations. 

Additionally, since the U(1)d effectively makes the dark halo a plasma (albeit a very cold, tenuous one), there may 
be other effects on structure formation that constrain this model [3H]. We have estimated that the timescale for the 
Weibel instability in our model is short compared to relevant timescales for galactic dynamics. If this instability has 
a dramatic effect when subhalos collide during the assembly of a galactic halo, our U(l)n could be excluded for the 
entire range of interesting parameters. Further work is required to before we reliably understand the quantitative 
effects of such instabilities on galactic dynamics. 

This work opens a window to new phcnomcnological possibilities within the dark sector. One avenue for further 
investigation would be the possibility of "dark atoms," which would arise if there were two different stable species 
with dark charge, each with an asymmetry in the number density of positive and negative charges (with one balancing 
the other to maintain overall charge neutrality). From there, one is free to contemplate dark chemistry and beyond. 
Dark matter constitutes a large majority of the matter density of the universe, and there is no reason to assume a 
priori that physics there is any less rich and interesting than that of ordinary matter. 
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